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1.0 INTRODUCTION 


Consider the differential equation 

r"(t) = G(r(t)) + u(t) (1-1) 

describing the motion of a spaceship, where r is the position vector, G(r) 
is the gravitational acceleration, and u denotes the thrust acceleration. 

Let us introduce a six-dimensional vector y = (r,v) where v = r'. 

Then equation (1-1) can be written in the form 


y' = f(t,y,u) 


( 1 - 2 ) 


where 


f(t,y,u) = (v, G(r) + u) 


(1-3) 


Since the velocity of any motion does not exceed the speed of light and 
the gravitational acceleration for motions (in a central gravitational field) 
which are bounded away from the center can be considered as bounded, the 
resulting function f can be assumed to be bounded, at least for physical 
motions, and the gravitational acceleration can be written in the form 


G(r) = -k(r)r 


(1-4) 


where 


k(r) = p/r^ if r > rg 


(1-5) 


and 


k(r) = p/rg^ if r < rg 


( 1 - 6 ) 


This modification of the gravitational acceleration will not change the mo- 
tion provided that the motion is outside the sphere with radius rg with 
its center at the origin. 
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It is easy to prove that after such a modification the function G 
is Lipschitzian and its smallest Lipschitz constant (as discussed by 
V. M. Bogdan and V. R. Bond in a paper entitled "A priori Global Estimates 
of Deviation of Solutions of Differential Equations Due to Perturbation With 
Applicatons to Analytical and Celestial Mechanics") is given by the formula 


I Gil = (2y)/ro3 


(1-7) 


Notice that in the above argument it is not essential that the central 
gravitational field has its center at the origin. The center of the field 
could be in motion and the above regularization procedure would still apply 
provided that the distance between the center and the spaceship is always 
larger than r^. Notice also that if the gravitation acceleration function 
is the sum of several central Newtonian fields then the above procedure 
applied to each of them will regularize the function G to make it 
Lipschitzian. For details see the paper referenced in the previous paragraph. 

Let us assume that the automatic control is given by the formula 


u(t) = g(t, y(t), tj, y(tj)) (1-8) 

where denotes the time that the telemetric reading of the position- 

velocity vector y(tj^) is taken. We assume that this information is trans- 
mitted to the control center, which is located at the origin of an inertial 
coordinate system. At the control center, this telemetric information is 
used to produce a computer program which is then transmitted to the space- 
ship. The program represents a function such that, when the program reaches 
the spaceship, new telemetric readings of the time and the vector y(t) are 
taken, and the program is executed, one obtains as a result the value of the 
control vector given by formula (1-8). We assume that this process is per- 
formed continuously, as by an analog computer. We shall be interested in 
such motions when the transmission time is not negligible, as when the space- 
ship is far away from the control center. 

Our first objective will be to find the time tj and the position- 
velocity vector y(t^) as an operator of the trajectory function y. 

The mathematical terminology and notions to be used in this paper can 
be found in the books Optimal Control by Athians and Falb (ref. 2) and 
Foundations of Modern Analysis by Dieudonne (ref. 4). For extensive lit- 
erature on delay differential equations and related topics, see the survey 
by Corduneanu and Lakshmikantham (ref. 3). For regularization of the gravi- 
tational acceleration used in this paper and explicit formulas for computing 
the smallest Lipschitzian constant for a vector function, see the aforemen- 
tioned paper by Bogdan and Bond. Valuable information on delay differential 
equations is contained in Driver (ref. 5) and Hale (ref. 6). Papers by 
Anderson and Bogdanowicz (ref. 1) and Lakshmikantham, Leela, and Moauro 
(ref. 7) are closely related to the results presented in this paper. 
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2.0 TIME DELAY OPERATOR 


Assume that a trajectory is defined by 


y(t) = (r(t), v(t)) 


(2-1) 


for all values of time. Let us fix time t and denote by T) = t - tj the 
delay caused by the transmission of the telemetric data to the control center 
and of the program back to the spaceship. This delay can be represented by 
the formula 

n = p + |r(t)l/c + )r(t - T))l/c (2-2) 

where p denotes the processing time required by the computers and c 
the velocity of light. (See the following sketch.) 



3 



Let us assume that for some number 


0 < q < 1 (2-3) 

every admissible physical motion is such that 

|v(t)| < cq (2-4) 

for any value of t. 

Under this assumption, we have the following estimate 
n < p + |r(t)|/c + |r(t) - v(s) ds|/c 

— I I I t-n ' 

< p + 2|r(t)|/c + nq (2-5) 

which yields the inequality 

n < (1 - q)“Mp + 2|r(t)|/c) = a(t) (2-6) 

Consider a function s defined by 

s(n) = P + |r(t)|/c + |r(t - n)|/c (2-7) 

for all n ^ 0* Observe that the function s maps the closed interval 

1 = <0,a(t)> (2-8) 

into I and it represents a contraction of this interval. Indeed, we have 
the inequality 

|s(n) - s(p)| < (1/c) |v(x)| dx < q|n - p[ (2-9) 

which proves that the function s is a contraction. Thus, there exists a 
unique number r\ satisfying equality (2-2) for every trajectory function y 
satisfying condition (2-4). 
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This means that the time tj can be treated as an operator of the trajectory 
function y and is given by the formula 


tl = TjCyKt) = t - n (2-10) 

3.0 OPERATORS ASSOCIATED WITH AUTOMATIC CONTROL OF SPACESHIPS 

In this section, we will investigate the operator and the operator 

defined by the formula 

T2(y)(t) = y(Ti(y)(t)) (3-1) 

for t belonging to the interval I and y belonging to a suitable space 
of continuous functions. Using these operators, we can convert the original 
spaceship control problem to the following form: 

y*(t) = f(t, y(t), u(t)) 

> for all t ^ 0 

u(t) = g(t, y(t), T]^(y)(t), T 2 (y)(t)) ) 

y(t) = x(t) for all t < 0 

where the function x is assumed to be known. 

Notice that a physical motion must satisfy the following Lipschitzian 
constraint 

|v(t) - vCt^)] < a|t - tj (3-5) 

since the acceleration of any spaceship cannot exceed a certain limit at which 
destructive stresses will be produced. Also, the fact that the velocity of 
the spaceship cannot exceed the velocity of light yields another constraint 

|r(t) - r(t]^)l < c|t - tj^l (3-6) 

for any t and 


(3-2) 

(3-3) 

(3-4) 
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In the sequel, we shall denote by I an interval of the form 


<0,b> or <0,oo) 


(3-7) 


Let Y be a Banach space and let x be a Lipschitzian function with 
Lipschitz constant w from (- 00 , 0> into Y. We shall denote by 


Lxw(Y) (3-8) 

the space of all functions y from the interval (- 00 ^ 00 ) into the Banach 
space Y such that 

|y(t) - y(t^)| < w|t - t^l (3-9) 

for all t,tj e I and 

y(t) = x(t) (3-10) 


for t < 0. 

Notice that this set of functions restricted to I can be considered as 
a subset of the space C(Y) of all continuous functions y from the interval 
I into the Banach space Y. In the space C(Y), we shall introduce a family 
of norms given by the formula 


y|^ = sup {e ^^|y(t)| I t G l} < «> 


(3-11) 


for any fixed real number k. Observe that from this definition follows the 
inequality 

ly(t)l < lylkel^t (3_12) 

for all tel. 

Notice that when the interval I is compact all the norms (3-11) are 
equivalent. When the interval I is unbounded, each value of k yields 
a different subset C^(Y) of the space C(Y) of all continuous functions. 
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Let Zj , Z 2 be two Banach spaces. Let us assume that D is a subset 
of the space Z^. We shall denote by 

L(D, Z 2 ) (3-13) 

the space of all operators U from the set D into the Banach space Z 2 
satisfying the Lipschitz condition; i.e., there exists a number L such 
that 


|u(z) - U(z]^)| < l|z - z^\ (3-14) 

for all z and in the set D. We shall denote by j|u|| the greatest 

lower bound of all the numbers L satisfying inequality (3-14). The value 
||u|| considered as a function of the operator U forms a seminorm on the 
space L(D, Z 2 ) * 

Lemma. Let Zp Z 25 Z 3 be some Banach spaces and let denote a 

subset of the space Z^ for i = 1, 2. If e ^(^1? ^i+ 1 ^ i = ly 2, 

then the composition operator U = i^sio^igs to the space L(Dp Z 3 ) and 

the following inequality 

IIU 2 U 1 II £ IIU 2 IIIIU 1 II (3-15) 

holds. The proof of the lemma is obvious. 


4.0 OPERATORS OF THE EXPONENTIAL TYPE 


Assume that given are two Banach spaces (i=l,2). For the sake of 

brevity, let us denote by = CCY^) the Banach space of all continuous 
functions from the interval I into the space Y^. Let denote a subset 

of the space for i=l, 2 . 


4.1 Definition 


An operator U will be said to be of the exponential type if and only 
if it belongs to some space L(DpC 2 ) of Lipschitzian operators and there 
exist a number b .and a positive constant such that ||u||^ < b for 

all k > k^ where || || ^ denotes the seminorm induced on the space L(DpC 2 ) 
by the norm | defined on the spaces (i=l, 2 ) of continuous functions. 
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In the sequel, we shall assume that we are dealing with motions y 
from the interval (-00,00) into the space Y = such that 


y(t) = (r(t), v(t)) 


(4-1) 


and 


y(t) = x(t) (4-2) 

for times t < 0, 

4.2 Theorem 


The operator T^^ considered from the set L^^(Y) into C(R) is of the 
exponential type for any continuous function x from (- 00 , 0> into Y having 
Lipschitz constant w < c. Moreover, 


((Tillk < 2/(c - w) (4-3) 

for all k ^ 0 . 

Proof. Take a fixed argument t g I and two functions 


Yi e Lxw(Y) (4-4) 

for i = 1, 2, and let 

Yi = (ri, v£) (4-5) 

be the decomposition into the radius vector component and the velocity 
vector. Let us introduce the notation 

Ti(Yi)(t) = t - rii (4-6) 

for a fixed value of t. Recall that the numbers rii satisfY the equations 

ni = (|ri(t)| + |ri(t - Tii)|)/c + p (4-7) 
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See equation (2-2) • 


Relations (4-6) and (4-7) yield the equality 

|Tx(yi)(t) - Tx(y2)(t)| = Im - 1 I 2 I (4-8) 

and thus we get the inequality 

c|ni “ TI 2 I 1 ' ^2^^ “ ^2^1 (4-9) 

Let us estimate the second term on the right-hand side of inequality 
(4-9). We have from the triangle inequality 

|ri(t - nj) - T 2 (t - ri2)I < - m) - rj(t - ri2)| 

+ |r^(t - ri2^ “ ^ 2 ^^ ~ ^^2^ (4-10) 

From the Lipschitz condition, the first term on the right-hand side of 
inequality (4-10) can be estimated by 

lr^(t - m) - r^Ct - ri2)| 1 w|ni - ri2l (4-11) 

Notice that the second term on the right-hand side of inequality (4-10) is 
equal to zero if 

t - ri2 1 0 (4-12) 

and it can be estimated in any case by 

|ri(t - ri2) - r2(t - 02 )! 1 ki - r2lke’^^^ '^2^ 

1 bl “ y2lk®'^'^ (4-13) 

as follows from the definition of the norm | |j^ (see formula (3-12)). 
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since we also have the estimate 


|ri(t) - r2(t)| < lyi - 72^6^*^ 


(4-14) 


for all t G I, we get from estimates (4-10), (4-11), and (4-13) and inequal- 
ity (4-9) the following estimate 


(c - w)|ni - TI 2 I 1 2lyi - y2lk®^*^ (4-15) 

for all t in the interval !• Multiplying both sides of inequality (4-15) 
by and dividing by c - w, we get 


e - TI 2 I 1 b|yi - 721k 


(4-16) 


for all t in the interval I, where b denotes the constant 


b = 2/(c - w) 


(4-17) 


Taking the supremum over all tel on the left-hand side of inequality 
(4-16), we get the inequality 

|Ti(yi) - Ti(y2)|k < b|yi - 72|k (4-18) 

for all y^, Y 2 from the set L^^(Y) and for k _> 0. Hence, we have the 
estimate 

l|Ti||k<b (4-19) 

for all k ^ 0. This completes the proof. 
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4. 3 Theorem 


The operator T 2 defined by formula (3-1) is of the exponential type 
from the set into C(Y) for every fixed function x with Lipschitz 

constant w < c. Moreover, 


ll^lllk < (c + w)/(c - w) (4-20) 

for all k > 0. 

Proof. Take any fixed tel, yj, J2 £ Ljj^,(Y), and notice the inequality 

|T2(yi)(t) - T2(y2)(t)| = |yi(Ti(yi)(t)) - y2(Tj(y2)(t)) I 

< |yi(Ti(yi)(t)) - yi(Ti(y2)(t)) I 

+ |yi(Ti(y2)(t)) - y2(Tj(y2)(t)) I (4-21) 

The first term on the right-hand side of inequality- (4-21) can be 
estimated from the Lipschitz condition as follows 


lyi(Ti(yi)(t)) - yi(Ti(y2)(t)) I < w|Ti(yi)(t) - Ti(y2)(t)| 

< w|iTlllk|yi - 72^6^^*^ 

< (2w/(c - w))|yi - Y2lk^^^ (4-22) 


The second term on the right-hand side of inequality (4-21) 


|yi(Ti(y2)(t)) - y2(Ti(y2)(t)) I 

is equal to zero if the value Tj(y 2 )(t) is negative and it can be estimated 
in any case by the quantity 

IyI “ y2lk®^*^ (4-23) 

as follows from the definition of the norm ] (see formula (3-12)). 
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From inequality (4-21), we get the estimate 


T2(yi) - T2(y2)lk i b|yi ” yzik 


(4-24) 


where 


b = (c + w)/(c - w) 


(4-25) 


Hence , 


||T2||k < b (4-26) 

for all k 0. This completes the proof. 

5.0 LIPSCHITZIAN FUNCTIONS INDUCE OPERATORS OF THE EXPONENTIAL TYPE 

Let us assume that Y, Z denote two Banach spaces. Let f be a con- 
tinuous function from the product I x Y into the space Z such that 

|f(t, yj) - f(t, y2)j 1 L|yi - y2| (5-1) 

for all tel and all Vi>y2 ^ Moreover, let us assume that the func- 
tion satisfies the inequality 


f ( t, y) I < L 


(5-2) 


for all tel and all y e Y. The infimum of all such constants L appear- 
ing in inequalities (5-1) and (5-2) will be denoted by ||f||. Denote by 


LipdxYjZ) 


(5-3) 


the set of all such functions. It is easy to prove that this set of func- 
tions with the norm ||f|| as defined above forms a Banach space. 
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For a positive number w, we define the set 


B„(Z) = {z e C(Z) ; |z(t)| < w} ( 5 _ 4 ) 

Lemma. If f e Lip(lxY,Z), then the operator F, defined by z = F(v) 
if and only if 

z(t) = f(t,y(t)) (5-5) 

for all tel and y e C(Y), is of the exponential type from the space 
C(Y) into the set B^(Z), where w = ||f||. Moreover, the estimate 

Mkl llfll (5-6) 

holds for all k ^ 0. 

Proof of the lemma. Notice that the operator F maps the set C(Y) 
into B^(Z). To prove that it is of the exponential type, take any two 
functions 

yi»yz e c(Y) ( 5 _ 7 ) 

From the inequality 

|f(t,yi(t)) - f(t,y2(t))| < ||f|| ly^Ct) - y2(t)| 

1 l|f|| Iyi - (5-8) 

valid for all t belonging to the interval I, we get the inequality 

|F(yi) ~ I k 1 II f II Iyi - Y2lk (5-9) 

for all functions yi>y2 £ C(Y) and all k ^ 0. This completes the proof. 
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6.0 PROPERTIES OF A VOLTERRA INTEGRAL OPERATOR 


We shall define the operator J by the condition 


z = J(y) 


( 6 - 1 ) 


if and only if 

z(t) = x(0) + /q y(s) ds (6-2) 

for every tel and 

z(t) = x(t) (6-3) 

for every t £ 0 and every y e 

Theorem. The Volterra integral operator J is of the exponential type 
fran every set B,,(Y) into the set Lxw(Y) for any fixed function x with 
Lipschitz constant w. Moreover, 


||j||k < 1/1^ 


(6-4) 


holds for all k > 0. 

Proof. The operator J is well defined from the set B^^(Y) into the 
set Lxw(Y). To prove that this operator is of the exponential type, take 
any two functions yi>y2 e B^,(Y). Notice the inequality 

|J(yi)(t) - J(y 2 )(t)| < /J lyi(s) - y 2 (s)| ds 

< /J lyi - Y2|ke'^® 

< (l/k)e*^*^ |yi - y2lk (6-5) 

Hence, we get the inequality 

|J(yi) - J(y2)|k < 1^1 " y2lk 

for any two functions yi>y2 frorn the set B^(Y). This completes the proof. 
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7.0 DYNAMICAL SYSTEMS WITH AUTOMATIC CONTROLS OF THE EXPONENTIAL TYPE 


Consider a dynamical system with the state equation 


y'(t) = f(t,y(t),u(t)) (7-1) 


where the control satisfies the equation 


u( t) = U(y) ( t) (7-2) 


for all tel. The following theorem establishes the existence and unique- 
ness of solutions of the system of equations (7-1) and (7-2) under the assump- 
tion that the function f is Lipschitzian and the operator U is of the 
exponential type. 

Theorem. If f e Lip(lxYxZjY) and the operator U is of the expon- 
ential type from the set Lj^^(Y) into C(Z) for every function x from 
(- 00 , 0> into Y with Lipschitz constant w = ||f||, then for every such func- 
tion X there exists a unique solution to the automatic control problem ex- 
pressed by equations (7-1) and (7-2) such that y(t) = x(t) for t < 0. 

Proof. First, notice that every solution of the system of equations (7-1) 
and (7-2), if it exists, has to belong to the set where w = ||f||. 

To prove the existence of the solution, notice that the system of equa- 
tions with the initial condition y(t) = x(t) for all t < 0 is equivalent 
to the following integral equation ~ 

y(t) = x(0) + f(s, y(s), U(y)(s)) ds (7-3) 

to be satisfied for all tel. To analyze the integral equation, let us 
introduce the operator X from the space Lj^^(Y) into the space C(Yxz) 
defined by means of the formula 


X(y) = (y, U(y)) (7-4) 


for all y e Ljj^^(Y). Notice that the operator X is of the exponential type 
on every set L^w(Y). 
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Using the Volterra operator J and the operator F induced by the 
function f over the space C(YxZ), we can write integral equation (7-3) in 
the equivalent form 


y = JFX(y) 


(7-5) 


Notice that the operator JFX maps every set L^^(Y) into itself, where 
w = ||f||> that it is of the exponential type. Moreover, we have the 

estimate 


llJFXlIk < (1/k) llfll b 


(7-6) 


for all k > 0, where b denotes a constant such that 


lx| 


k 


< b 


(7-7) 


for all k > kg. Thus, if we select a value of k such that 

k > max{ b|[f(|, kg } (7-8) 

the operator JFX will become a contraction map as follows from formula (7-6). 
Since the set L^„(Y) can be considered as closed in the Banach space Ck(Y), 
it represents a complete metric space. Hence, by the Banach contraction 
mapping theorem there exists a unique fixed point y £ Mjj^^(Y) for the 
operator; i.e.. 


y = JFX(y) 


(7-5) 


This completes the proof of the theorem. 


8.0 APPLICATION TO AUTCMATIC CONTROL OF SPACESHIPS 


Consider the original problem. Let us assume that we are given the 
motion of a spaceship 


y'(t) = f(t,y(t),u(t)) = (v(t), G(r(t)) + u(t)) 


(8-1) 
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x^dth controls given by 


u(t) = g(t, y(t), TjCyKt), T 2 (y)(t)) (8-2) 

We assume that Y = and Z = R^. The operators Tj^,T 2 are defined as in 

the previous sections. We assume that the motion of the spaceship y(t) = x(t) 
is known for t < 0. 

Theorem. If f e Lip(lxYxZ,Y) and g e Lip(lxYxRxY,Z) and the oper- 
ators Tj and T 2 are defined as before, then for every function x from 
(- 00 , 0> into Y with Lipschitz constant w = ||f|| < c there exists a 
unique solution to the automatic control problem given by equations (8-1) 
and (8-2) such that 


y(t) = x(t) 


(8-3) 


for t _< 0. 

Proof. The operator U, defined by the condition 


u = U(y) 


(8-4) 


if and only if 


u(t) = g(t, y(t), T^(y)(t), T 2 (y)(t)) (8-5) 

for all tel, maps the set L^^(Y) into the space C(Y) of continuous 
functions and is of the exponential type. To prove this, take any two func- 
tions yi,y2 ^ assume that U£ = UCy^) for i = 1, 2. Then we 

have the following estimates 


|ui(t) - U2(t)| < ||g|| (|yj(t) - y2(t)| + |T^(yj)(t) - T^(y2)(t)| 

+ |T2(yi)(t) - T2(y2)(t)|) 


< blyi - y2lks'^‘^ 


( 8 - 6 ) 


for all tel and yi»Y2 ^ ^xw^^^> where 

b = ||g|| (1 + 2/(c - w) + (c + w)/(c - w)) 


(8-7) 
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Multiplying the expressions in inequality (8-6) by e and taking the 
supreraum over all t e we get the inequality 


I’^i ~ “2lk 1 ^\y\ - yi\k (8-8) 

for all k > 0 and yijY2 G The last inequality proves that 

Ihllk 1 b (8-9) 

for all k > 0; i.e., the operator U is of the exponential type from every 
set L^^(Y) into the space C(Z) of continuous functions. 

The above proves that the conditions of the theorem in section 7.0 are 
satisfied, and thus for every initial condition 


y(t) = x(t) for t < 0 


(8-3) 


there exists a unique solution to the automatic control problem given by 
formulas (8-1) and (8-2). This concludes the proof of the theorem. 

9.0 CONCLUSION 


The main result of this paper is the proof of the existence and unique- 
ness of the solution to the automatic control problem with nonlinear state 
equation of the form 


y* = f(t,y,u) (9-1) 

and nonlinear operator controls 

u = U(y) (9-2) 

acting onto the state function y which satisfies the initial condition 

y(t) = x(t) (9-3) 

for t < 0. 
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The required regularity conditions on the function f and the operator 
U are the following: The function is of Lipschitzian type in the variables 

y and u and is bounded and continuous on its domain. The operator U is 
Lipschitzian from every set into the space of continuous functions 

for every initial function x with Lipschitz constant w = ||f||. Moreover, 
the induced seminorm ||u||^ as a function of k should be bounded for suf- 
ficiently large k. This class of nonlinear operators is said to be of the 
exponential type. Generalization of these results will appear in reference 8. 
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